Collisional cooling of ultra-cold atom ensembles using Feshbach resonances 
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We propose a new type of cooling mechanism for ultra-cold fermionic atom ensembles, which 
capitalizes on the energy dependence of inelastic collisions in the presence of a Feshbach resonance. 
We first discuss the case of a single magnetic resonance, and find that the final temperature and 
the cooling rate is limited by the width of the resonance. A concrete example, based on a p-wave 
resonance of ^''K, is given. We then improve upon this setup by using both a very sharp optical 
or radio-frequency induced resonance and a very broad magnetic resonance and show that one can 
improve upon temperatures reached with current technologies. 



The technology of cooling atomic ensembles has been 
one of the most important developments in physics over 
the last decades Cooling atomic samples has two in- 
gredients. The first is a "knife" that selectively removes 
atoms with the largest kinetic energy. Secondly, elastic 
collisions between the atoms thermalize the remaining 
atoms. It has been a critical ingredient in creating Bose- 
Einstcin condensates in improving atomic clocks 
and studying atomic properties [3]. The temperatures 
that have been achieved in bosonic gases are now well 
below a nano Kelvin. For fermionic systems, however, 
temperatures are just below micro Kelvin, due to limita- 
tions of current technology. There is therefore interest in 
developing improved cooling methods for fermions. 

We propose a new type of cooling mechanism that uses 
inelastic scattering processes due to a narrow magnetic 
Feshbach resonance 0, [3|- Here, a molecular state that 
is resonantly coupled to two scattering atoms acts as a 
"knife" that changes the internal state of colliding atoms 
in a strongly energy-selective manner. These states are 
then either untrapped or have gained sufficient kinetic 
energy to quickly leave the trap. Elastic scattering near 
the resonance leads to thermalization. For a resonance 
to generate loss, the internal state of the atoms cannot 
be the lowest energy state, as there need to exist states 
into which they can scatter. For samples of atoms in their 
lowest state, we further propose combining a narrow opti- 
cal [§] or radio-frequency (rf) induced resonance with 
a broad magnetic resonance. Here, the narrow resonance 
generates the loss processes, whereas the broad magnetic 
resonance drives thermalization. Wc will show that we 
can reach lower temperatures than feasible with a single 
resonance. 

Figure [U shows a schematic representation of the cool- 
ing process. We assume a gas of Fermi atoms in a single 
state, although the ideas can be generalized to bosonic or 
multi-species and -states fermionic gases. In the degen- 
erate regime the atoms with mass m form a Fermi sea, as 
shown in the figure, with Fermi energy Ep = kp/{2m), 
Fermi momentum kp, and a temperature T less than 
Ep/ks- We consider a narrow Feshbach resonance, that 
for two colliding atoms with momenta k and p induces 
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FIG. 1: Schematic representation of the loss processes that 
lead to cooling. We show the momentum distribution {k^, ky) 
of a degenerate Fermi gas (blue, inner feature), and the in- 
elastic scattering rate Kin{k,p) (orange, ring-shaped feature), 
both in arbitrary units, with a scattering partner with a mo- 
mentum slightly above the Fermi surface at p = {px,Py) = 
(1.05fcF, 0), and with a resonance energy that corresponds to 
a relative momentum of l.Obkp- The losses are largest where 
the two surfaces approach each other. This corresponds to 
atoms with momenta located on opposite sides of the Fermi 
sea. (The momenta along z are not shown for clarity.) 
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This coefficient is only a function of the relative collision 
energy E and is strongly localized around the resonance 
energy iS^es, which can be controlled by external mag- 
netic field. The energy E is given hy E = fc^/(2m,,) = 
where k^ = (fc — p)/2 is the relative momen- 
tum and rrLr = m/2 is the reduced mass. Finally, the 
total energy width is Ttot = Fq -I- r(i?), where Tq/H is 
the linewidth of the resonant state and T{E)/h is the 
collision-energy-dependent stimulated width. 

In order to generate cooling, we need to choose the res- 
onance energy so that only atoms with momenta larger 
than kp are lost. The largest relative momentum is about 
kp, corresponding to a pair of atoms on opposite sides of 
the Fermi sea. Consequently, Eres must be set slightly 
above kp/{2mr) = 2Ep. The exact amount that the 
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resonance energy lies above the Fermi energy will be de- 
termined by the temperature of the gas and Ttot- More- 
over, the resonance energy needs to be gradually lowered 
in time by changing the magnetic field, as the atom num- 
ber decreases due to the losses, and thus the Fermi energy 
decreases. 

At the beginning of the cooling process the tempera- 
ture is much larger than Ttot/ks- We will find that the 
smallest temperature that can be achieved is 4:Ttot/kB- 
This suggests the use of arbitrarily narrow resonances. 
However, the thermalization rate is also influenced by 
the resonance. In fact, the total elastic scattering rate 
coefficient Kei{k,p) is 

Typically, during the cooling process it is preferrable that 
the Fermi gas is close to thermal equilibrium, and there- 
fore we require the ratio K^j/Kin = r(£')/ro S> 1 for 
E « 2Ef- Consequently we have Vtot ~ ^{E). Since we 
want a temperature that is as low as possible, this leads 
to competing requirements on Fq and T{E). 

For fermionic atoms in the same internal state, only 
odd partial wave scattering exists. In fact, for ultra-cold 
atoms we only need to include the poi £ ^ 1 partial wave. 
Moreover, the energy dependence of F(£') is enforced by 
the Wigncr threshold laws. Here, this leads to F(£') = 
AE^^^, where A is an intrinsic property of the resonance. 

In Ref . , a p-wave resonance was characterized in the 
collisions of fermionic atoms in the hyperfine state 
|9/2,— 7/2). For the mi = component of the p-wave 
resonance, they find Fo/fcs = 0.9 nK and T{E)/kB = 
1.4 X 10-3 nK at E/ks = 1 nK. 

Classical limit. Before we consider cooling a Fermi gas 
in the degenerate regime, we treat the classical limit of 
a three-dimensional homogeneous gas, in which we as- 
sume that the momentum distribution /(p, t) remains a 
Maxwell-Boltzmann distribution throughout the cooling 
process [ll]. In fact 

= tJ7(n,T)cxp(V/(2™fcBr)), (3) 

where w{n,T) = {2Trh)^n/{2TrmkBT)^^'^ is the phase 
space density, and only the particle density n and the 
temperature T are time dependent. Its time evolution is 

/d^k 
^—^K,nik,p)f{k,t)f{p,t). (4) 

In the hmit that Fq + r{E) < ksT, we find that the 
approximate time evolution for n and T is given by 

dtn = -7,;„(t)n, dtT = -j,n{t) ("^^ - ^) (5) 
with the rate 

7,„W - ^23/^n7(n,T) Jf{.f-\ e-^'-/(^-^),(6) 



being linear in n, and time dependent through n, T and 
Eres- For this approach to be consistent, the thermal- 
ization rate during the evolution needs to be larger than 
the rate ^in{t). Within the classical theory we find that 
the thermalization rate is equal to Eq. [S]with Fg replaced 
by T{E) in the numerator. Therefore we have to require 
T{E)/Tq ^ 1, and 7m(i) becomes independent of F(i?). 

As an example for the dynamics described by Eqs. ([5|). 
we consider a process in which the resonance energy 
tracks the temperature at a fixed ratio A, i.e. Eres{t) = 
\kBT{t). We find the solution 

n{t) = no(l-<Ad)^, T{t) = To{l-t/t,,i)i^ (7) 

where hq and Tq are the initial density and temperature, 
and the classical cooling time td is given by 

ox 7 

l/tci - -^-/Ut = 0). (8) 

The phase space density increases as w{n,T) = wo{l — 
t/tci)^^, where zuq is the initial phase space density. 
When this approaches one, the system reaches degen- 
eracy and the classical limit breaks down. This occurs 
shortly before td, if wq ^ 1. 

For the ^'^K example, A = 11/2, Tq = 1 /^K, and 
uq = lO^'^ cm"'^, we find td « 2 s. This is a realistic 
time scale for current experimental setups, motivating 
the subsequent, more in-depth analysis. 

Single resonance cooling. We now consider cooling of 
a gas of spin-polarized fermions in the degenerate regime 
using a single Fcshbach resonance. We use the quantum 
kinetic theory of Rcfs. [l^ . The quantum dynamics 
is then fully given by the evolution of the momentum 
state occupation, which satisfies a homogeneous quantum 
Boltzmann equation. We further assume a spherically 
symmetric momentum distribution, and the momentum 
distribution / becomes a function of kinetic energy e = 

/{2m) only. This fimction /(e) satisfies 

p{ei)dtf{ei) = -K(ei)p(ei)/(ei)-/(ei) (9) 

where the density of states per unit volume is p{e) = 
47rm3/^'\/2e/(27r?i)'^, the loss rate k is given by 

A^(ei) = P(e2)de2/(e2) j dcoseK,„{E). (10) 

We denote Kin{E) = Kin{ki, k2), as the inelastic loss 
rate coefficient only depends on the relative energy E = 
ei/2 -I- 62/2 — Y^eleJ cos6', and 9 is the angle between ki 
and k2. 

The collision integral is given by 

/(ei) J ^62^631^(61,62,63,64) (11) 

{/(ei)/(62)(l-/(63))(l-/(64)) 

-(l-/(ei))(l-/(62))/(63)/(e4)} 
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FIG. 2: CoUisional cooling of a degenerate Fermi gas with 
a Feshbach resonance as a function of time for the example 
described in the text. Panel a) shows the distribution /(e) 
as a function of time. Panel b) shows the temperature T and 
ATtoti^Ep)- Finally, panel c) shows T/Tp and the density n. 



where 64 = ei + 62 — 63, and the coUision kernel W is 
^61,62,63,64) = -J^^ _/p dPa{E). (12) 

The clastic cross-section (j{E) = Kei{E)/vr{E) and we 
use Kf,i{E) = Ki,i{ki,k2) as it only depends on the 
relative energy E = ei + 62 — P^/(4m), where P is 
the total momentum in the collision. The integration 
bounds are given by Pmm = max(|pi -P2I, \ Pz--P 4\) and 
P^nax = min(pi -l-p2,P3 +P4), where pi = ^/Sroel. 

We can estimate the thcrmalization rate of the system 
in the quantum degenerate regime, based on the Boltz- 
mann equation 0. Following the procedure outlined in 
[ist , we linearize the Boltzmann equation around a Fermi 
distribution /o(e) with Fermi energy Ep and temperature 
T, that is, /(e) ^ /o(e) + /o(e)(l - /o(e))V'(e), where 
i^{e) is a small deviation and the functional form ensures 
that the fluctuations are localized around the Fermi en- 
ergy. Then we find that the thcrmalization rate for states 



close to the Fermi energy is given by 

— ^ ^^f^W{EF.EF.Ep,Ep), (13) 

refiecting that the only contributions to the collision in- 
tegral of Eq. [11] are processes close to the Fermi energy. 

The value of W{Ep,Ep, Ep, Ep) can be estimated by 
realizing that the integral in Eq. [T^ runs from zero to 
2kp, and therefore the relative energy from zero to 2Ep. 
As the resonance energy in our cooling scheme will be 
larger than 2Ep by an amount of the order of r(2Ep), 
we find that W oc yJV{2Ep) / Ep and the thcrmalization 
rate is 

— ^ [keTfA^/^E-'^^lh. (14) 

Tth 

The quadratic temperature dependence is typical for a 
Fermi gas in the degenerate limit. For the losses we sim- 
ilarly find a time scale l/n ^ TaT{2Ep)/{hEp). 

We now solve the quantum Boltzmann equation nu- 
merically to study our cooling process, starting from a 
Fermi distribution for the atoms. The initial resonance 
energy E^es is set well above twice the Fermi energy. We 
then gradually lower E^.^^ to eliminate atoms with large 
kinetic energy. 

The final i?res, and thus Ep, and the time scales of the 
cooling process can be estimated from our expectation 
that the smallest T/Tp is of the order Vtot/Ep. Mini- 
mizing this with respect to Ep gives Ep = {2Tq/ A)^/^ . 
The time scales Tth and t; will be largest at that final 
value. They are approximately Tth ^ ^/(^^g^^) and 
Ti^h/{AV^Ty'). 

Figure [2] shows an example of cooling for the "^^K 
resonance described before. The atomic ensemble has 
an initial temperature of 1.0 /xK and a chemical poten- 
tial of fj,(0)/kB = 0.5 fiK, corresponding to an initial 
density of n w 5 • 10^"^ cm"-^. We choose Eres{0) ~ 
Eres.t 8 fiK, well above 2Ep « 2.5 fiK and T. We 
choose Eres.f = 0.2 /iK as the final value for the reso- 
nance energy, based on the above estimate for the op- 
timal T/Tp. We use an exponential time dependence 

Eres{t) {Eres,i ~ -E^es,/) exp(-i/io) + Eres,f, with a 

timescale <o = 5 s, that is larger than the estimates for 
Tth and T;. In Fig. [^^ we show the distribution /(e) as 
a function of energy and time. It becomes visibly colder 
in the process. Throughout the simulation /(e) is fairly 
close to thermal, as expected for r(2£'F)/ro ^ 1. Con- 
sequently, one can fit /(e) to a Fermi distribution, and 
assign a temperature and a chemical potential at any 
point in time. Figure [^b shows the fitted temperature as 
a function of time. It gradually approaches the expected 
temperature ^Ttot{2Ep)/kB- The fastest cooling rate or 
slope, around t « 3 s, is consistent with the rate ex- 
pected in the classical limit, Eq. [H The cooling is slower 
for lower temperatures consistent with Eq. 1141 The tem- 
perature has decreased by two orders of magnitude in 
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FIG. 3: CoUisional cooling with an optical resonance as a 
"knife" and a magnetic resonance for thermalization. We use 
themassof '"'K,r°P*(£;)//fcs = 40-10"^ nK[S/(l nK)]^/^^ and 
Y'o'^ /ks ~ 1.8 nK. For the time dependence of the resonance 
energy we use E°^l i = 4.1 fiK, E°ll f = 1-2 /iK, and to = 6 s. 
a) /(e,f). b) T/Tp'. 



the process, from 1 /iK to 10 nK. In Fig. [2j: we show 
T/Tp and n. T/Tp is reduced by a factor of 6, from 0.8 
to 0.15, while simultaneously the density has decreased 
from 5 X lO^^ ^ q.T x lO^^ cm^^ 

We have performed multiple simulations with various 
initial densities, temperatures, resonant energies and tg. 
We observe the same qualitative behavior, except for very 
small values of to, where we lose atoms too quickly and 
the system does not equilibrate. Most importantly, we 
find that fc^T approaches 4rtot. By optimizing the func- 
tional form of Eresit), this temperature could be achieved 
in a shorter time. However, we do not expect to reach 
significantly lower temperatures. 

Cooling with two resonances. The elastic and inelas- 
tic scattering rate coefficients are governed by T{E), 
which leads to contradictory requirements. To over- 
come this limitation, we propose the use of two p-wave 
resonances: a narrow optical or rf induced one, which 
acts as the 'knife', and a broad loss-less magnetic one 
which thcrmalizes. We locate the magnetic resonance 
such that the elastic rate coefficient is unitarity limited 
at Kei = VrTrh^ /k"^. To ensure a loss-less magnetic reso- 
nance the atoms must be in the lowest hyperfine state and 
the field driving the narrow transition must not couple 
to the Feshbach molecular state [31 • On the other hand, 
the loss from and the width of the optical (rf) resonance 
can be controlled by the laser intensity (rf field). This 
creates a very narrow "knife" with a negligible contribu- 



tion to the elastic rate coefficient. Its resonance location 
E°ll is lowered in time. 

We have repeated the analysis of the previous section. 
The thermalization rate is driven by the magnetic res- 
onance and given by l/rth ^ [ksTY /{hEp). Hence, 
unlike for the single resonance case, thermalization does 
not require T°p\E)/T°^'* > 1. In fact we find that 
V°'p\2Ef) k. r°P* leads to the fastest cooling for a given 
r°Q*. Figure [3] shows an example of the cooling pro- 
cess. The initial density is 4 x 10^"^ cm~'^, the initial 
temperature is 0.3 /iK. The system is cooled down to a 
final density 1.9 x 10^'^ cm~'^ and a final temperature 
0.017 /xK. This demonstrates that T/Tp « 2.5 x 10"^ 
can be achieved with this cooling process. 

In conclusion, we have proposed a new cooling mecha- 
nism which uses the energy selectivity of a Feshbach reso- 
nance. We first discussed the limit of a classical thermal 
gas, before we turned to a quantum kinetic simulation 
of a degenerate Fermi gas. The case of a single reso- 
nance shows cooling to the regime around T/Tp ~ 0.1 
for ^''K, for an appropriately chosen resonance. We then 
improve on this setup by using one narrow resonance for 
the loss process and a broad magnetic resonance for ther- 
malization. This setup can create temperature regimes 
competitive with current technology. 
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